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When chiral symmetry is spontaneously broken, the low-energy part of the Dirac 
operator spectrum can be computed analytically in the chiral limit. The tool is 
effective field theory or, equivalently in this case. Random Matrix Theory. 



1 Introduction 

Contrary to what one might naively have expected, it is possible to predict 
the analytical behavior of the Dirac operator spectrum in QCD and QCD-like 
theories for eigenvalues very close to zeroj^ This is the low-energy regime of the 
theory, and ordinary QCD perturbation theory is completely irrelevant here. 
In fact, the analytical predictions turn out to be not only non-perturbative in 
the QCD coupling constant, but exact. These developments began eight years 
ago with a paper by Leutwyler and Smilgau, where for the first time it was 
shown that if QCD undergoes spontaneous chiral symmetry breaking, there 
are constraints on the Dirac operator eigenvalues near the origin. So the idea 
is not new, but there have been many new developments in just the past few 
months. Here I summarize some of the recent work I have been involvecLin 
myself. For a comprehensive review of the subject up to last year, see ref. □. 

What originally pushed this subject forward was the suggestion by Shuryak 
and Verbaarschot H that the lowest part of the Dirac operator spectrum could 
be computed analytically by means of Random Matrix Theory. Much work 
has gone into understanding this remarkable connection between chiral sym- 
metry breaking in gauge theories and Random Matrix Theory. In particular, 
it has been shown that the chiral (flavor) symmetry breaking in gauge theories 
coupled to fermions can be classified pr£cjajely according to the main categories 
of Random Matrix Theory ensembles □'ErQ. In the limit relevant for applica- 
tions to the Dirac operator spectrunLthese Random Matrix Theory ensembles 
actually provide universality classesu: results do not depend in detail on the 
particular distributions of the random matrices. 

Random Matrix Theory should not be required to derive results for the 
Dirac operator spectrum; all should follow from field theory alone. It is one 
of the more recent discoveries that this is indeed the case. There are now 
three independent means of deriving analytical results for the Dirac operator 
spectrum: 
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• Random Matrix Theory 

• The supersymmetric method 

• The rephca method 

Both the so-cahed supersymmetrici and repHcaiB methods stand out, as 
they rely only on the field theory formulation itself. Because results are claimed 
to be exact, all three methods should yield identical results. This is indeed the 
case. One pleasing consequence is that one can now use any of the three 
different formulations according to what is most convenient. Generically the 
Random Matrix Theory formulation is the easiest for analytical calculations, 
but there are cases-where it is actualht simpler to use the effective field theory 
partition function liiJ (see also refs. E3't3). Even in those cases, the explicit 
derivations have usually gone through Random Matrix Theory at intermediate 
steps. 

In this short review I will cover three different topics on which there has 
been progress in just the last few months. These are: (i) the analytical com- 
putation nLindividual eigenvalue distributions til (and subsequent lattice mea- 
surementsllB), (ii) the replioa-tnethod as applied to the finite- volume effective 
partition function of QCD and (iii) lattice Monte Carlo measurements 
of the low-energy end of the Dirac operator spectrum near the chiral phase 
transition EZrEa. 

2 Distributions of individual Dirac operator eigenvalues 

Analytical expressions for the distribution of the first non-zero Dirac operator 
eigenvalue have been known for some time in the case of all three chiral ensem- 
bles. Recently a much more general analytcial expression was found for the 
probability distribution of the fcth sniallest Dirac operator eigenvalue (again 
counted above the exact zero modes) EJ. The formula is completely general, 
holds for any number of (possibly massive) fermions, and for any of the three 
major chiral universality classes, as labelled by their Dyson index (3. The only 
technical restriction is that for the universality class of /3 = 1 (corresponding 
to gauge group SU(2) and fermions in the fundamental representation), the 
formula only works for a sector of odd topological charge ly. While the analyt- 
ical expression has been derived in the framework of Random Matrix Theory, 
the final formula turns out to depend only on the effective field theory parti- 
tion functions. It generalizes a related formula for just the lowest eigenvalue 
distribution in the case of the (3 = 2 universality class, which, also expressed 
the distribution in terms of finite volume partition functions EJ. 
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Before presenting the resulting closed analytical expression for the distri- 
bution of the fcth smallest Dirac operator eigenvalue, let us first recall some 
basic facts about the connection to Random Matrix Theory. Because chi- 
ral symmetry is presumed spontaneously broken, the large-volume euclidean 
partition function is, in the chiral limit, entirely dominated by the pseudo- 
Goldstone bosons. The dominance is very strong: contributions from all other 
physical excitations are exponentially suppressed in the masses, with a coef- 
ficient in the exponent that grows linearly with the size of the box. As the 
size of the box is sent to infinity, the Goldstone bosons dominate the euclidean 
partition function entirely. It is in this sense that the results for the Dirac op- 
erator eigenvalues are exact: they can be made to reach any desired accuracy 
by taking the volume large enough, and the quark masses small enough. By 
imposing on the four-volume that V ^ 'm^'^ only the zero mode piece of the 
Goldstone field tt(x) survives, and the effective partition function, in a sector 
of topological charge v becomes lil 



Ju{Nf) [2 



(1) 



This partition function, a zero-dimensional group integral has, surprisingly, a 
Random Matrix Theory representation^. To get it, one starts withB 



Zi''^\{m,}) = 



Nf 

dW Yl (let {iM + mf) exp trV{M'^) 



where 



M 





W 



(2) 



(3) 



The complex matrix W in eq. ^ is of size (N + v) x N, and the potential is 
essentially not constrained beyond yielding a non- vanishing spectral density at 
the origin, p{0) ^ 00. The limit — > oo is taken in the Random Matrix Theory 
partition function in such a way that the combinations fii = miir p(Q)2N , are 

kept fixed. In this limit the partitjon functions z'^^^ and z\j^^^ become equal 
up to a /i^-independent constant^. 

For practical computations an eigenvalue representation of the above ma- 
trix integral is convenient. Taking a general Dyson index /3 = 1, 2 or 4, it can, 
up to an overall and irrelevant normalization factor, be written as 



N, 



N 



(4) 
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where Xi are the eigenvalues of Af^. Here the potential has been chosen to be 
just Vix) = X. This is precisely permitted because of universality of the final 
results 13. Due to symmetry under ;/ ^ — it is convenient to restrict oneself 
tov>0. 

The unnormalized joint probability distribution for all N eigenvalues is 



AT _ Nf N 

\Xi — Xj \^ 

1=1 j=l i>j 

(5) 

Similarly, the unnormalized joint probability distribution of the k smallest 
eigenvalues {0 < xi < • • • < Xk-i < Xk} is given by 

poo 

^7V,fc(2^i'---'^fe;{™?}) = / dxk+i---dxN PN\^i,---^XN;{m^}) 

k Nf k 



I Xi Xj I ^ 



i=l j=l i>j 

poo N Nf k 

X / dXk+l---dXN n {^^^''^'^^'''^'^^'I[i^^+^')]li^^-^ly 
•'^i' i=k+l j = l j = l 

N 

X n \^^-^Jf^ (6) 
i>j>k+l 

Shifting Xi ^ Xi + Xk in the integrand gives 

k Nf 

i=i j=i 

k i-oo N Nf 

x[]|x,-x,f / n (^dx,e-^^^x^{x,+Xkf(''+^'>/^-^l[{x,+m^^+Xk) 

k-1 N 

j=l i>j>k+l 

To finally get the probability distributions of the Dirac operator eigenvalues 
one takes the microscopic limit N oo with, in this particular convention, 
= TTp{0)y/xi = \/ 8N Xi and pj — 'Kp{0)mj = ^/SNrnj kept fixed. In this 
large- TV limit the difference between partition functions based on N — k and 
N eigenvalues disappears. One immediately sees that the new terms in the 



integrand of can be interpreted as arising from new additional fermion 
species, with the partition function now being evaluated in a fixed topological 
sector of effective charge v = 1 + 2/(3. The only restriction is that for /3 = 1 the 
topological charge v must be odd (because otherwise the number of additional 
fermion species will be fractional). 

Taking into account the definition (||), this gives: 



k 



1=1 

3=1 i=l 
Nf k-1 Nf ^iP) 



j = l i>j j"=l {it^if) 



In the partition function in the numerator each of the Nf fermion masses have 
been shifted according to A*? + Cfc- There are /3(fc — 1) new masses Ck~Ch 

i — 1, . . . , k — 1, each of them being (3 times degenerate. Finally there are also 
(3{i' + l)/2 — 1 new fermions, all of mass Cfc- The overall normalization factor 
C is fixed by the requirement of probability conservation. 

To get the individual eigenvalue distribution of the fcth eigenvalue, one 
simply integrates out the previous k — 1 smaller eigenvalues, viz., 



Jo J(l ■'Qk-2 



c^r(Ci,---,a-i,C;{Mj) ■ (9) 



k-l^k 

The general formulas (||) and may look rather complicated, but they 
actually simplify considerably in a number of interesting situations. For in- 
stance, in the physically, most interesting case of QCD (which belongs to the 
[3 — 2 universality classa), we get in a sector of topological charge v — Q: 

fc-1 Nf k-l 

c.,(Ci,...,Cfc;{M) = Ce-^l/\uX{(c^X{{C.l+^^,))X{{C^-CW■^ 

(10) 



^2 ({x^FTc^} , Va^„ VCl^^, . . . - Cfc^i, ^Cl CLi 



MM) 

The finite-vpljume partition functions involved here are known in closed ana- 
lytical formlla, 

Z,{{^l.,}) = detA{{^l})/A{p') , (11) 



where the determinant is taken over the Nf x Nf matrix 

MM) = ^Ji^^h+,-l{^i^) , (12) 

and 

Nf 

Mp') = Uif'l-t^-)- (13) 

i>j 

With this convention the normahzation factor is C = 1/2 for aU values of k, 
Nf and ly. 

The analyticaLformula ( p^ ) has very recently been tested by lattice gauge 
theory simulationstB, with quite remarkable agreement even at relatively small 
volumes. These simulations were done with staggered fermions, which are al- 
most totally insensitive to gauge field topology at the couplings we are con- 
cerned with hereE2l. This means that comparisons should be done only with 
the ly — analytical predictions. Shown in figure 1 is the result of a large- 
statistics computation in quenched SU(3) gauge theory. The individual Dirac 
operator eigenvalues indeed do have distributions that fall right on top of the 
analytical predictions. 
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Simulations with dynamical fcrmions in SU(3) gauge theory display the 
same degree_af accuracy (although for obvious reasons the statistics here is 
much lower) E3. Recall that both masses and eigenvalues are rescaled with 
the same factor oi I^V. This means that to get as good accuracy as possible 
one needs to probe the theory of nearly massless quarks. In figure 2 is shown 
an analogous plot for the theory with Nf = 1 (strong coupling) staggered 
fermions (at a mass value of m = 0.003 and a volume of 6^). 
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Fig. 2 Individual Dirac operator eigenvalue_distributions from lattice Monte Carlo simula- 
tions of SU(3) gauge theory with Nf = lllj, compared with the analytical predictions of 
eq. (0). 

One sees the same kind of detailed agreement as in the quenched theory. 
This establishes quite clearly that not only can the microscopic spectral density 
of the Dirac operator be computed analytically in the relevant scaling window, 
also individual Dirac operator eigenvalues are falling right on the analytical 
predictions. 

3 The replica method 

In lattice gauge theory simulations one often starts out with quenched calcu- 
lations: the fermion determinant is entirely ignored in the averages. This is 
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just a simple (sometimes accurate, sometimes not) approximation, introduced 
only in order to save computer time. A more sophisticated approach is "partial 
quenching" . Here one does simulate the full theory with dynamical fermions, 
but in addition one extracts more information from the same configurations. 
This is done by computing quenched averages in the field configurations that 
already included the effects of dynamical fermions. At first sight it may sound 
absurd to do this, but in fact there is much genuine physics to be extracted 
from such partially quenched simulations. In addition, it turns out that at 
the analytical level partial quenching is just what is needed in order to derive 
properties of the Dirac operator spectrum. In this case partial quenching is 
not at all considered as an approximation to the full theory; rather, it is used 
to obtain genuine physical information about the real theory, the one with 
propagating dynamical quarks. 

To make a partially quenched average, we must measure some correlation 
function of new fictituous "valence" quarks that do not affect the configu- 
rations. The perhaps most simple quantity is the partially quenched chiral 
condensate. We can get this condensate by adding valence quarks, differ- 
entiate with respect to their mass, and subsequently taking the limit Ny — > 0: 

= hm — - — InZ^ 



1 d d_ 



(14) 



In the last line we have formally expanded the partition function as a Taylor 
series in Ny. It is not obvious that precise meaning can be given to such a 
notion, since the function involved a priori is known only at integer values of 

Ny. 

This particular way of deriving the partially quenched chiral condensate 
is known as the replica methodj. Its applicability in certain condensed matter 
contexts has been questioned EJ, but last year there was considerable progress 
ixjwards understanding how to apply this method to derive spectral properties 
c3. In the present context of QCD, the partially quenched chiral condensate in 
t.Vip same finite- volume scaling regime as discussed above was computed in refs. 
Elilj. The replica method trivially works in ordinary QCD perturbation theory 
(it simply kills all closed fermion loops, as required). It is a more non-trivial 
fact that the replica method also is suited for deriving quantities that are non- 
perturbative in the QCD coupling constant. For instance, in |tlie low-energy 
framework of effective chiral Lagrangians the replica methpdc3 works quite 
analogous to the previously known supersymmetric methodB^. 
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Although in principle the partially quenched chiral condensate is an un- 
physical quantity, it turns out that this is not quite so. There is one tiny bit 
of this quantity that coritains important physics: this the discontinuity across 
the imaginary mass axis I, which gives the spectral density. In the microscopic 
scaling region: 

pt\c,{l^^}) = ^Disc 

= ^[^A^C + e; {f,,}) - E(*C - e; . (15) 

ZTT 

In ref. i two expansions of the partially quenched chiral condensate were con- 
sidered: small-mass and large-mass expansions. It turned out that neither 
were suitable for deriving the spectral density of the Dirac operator. The 
small- mass expansion suffered from so-called de Wit-'t Hooft poles Ea, while 
the large-mass expansion of ref. based as it were on the leading saddle 
point, only gave the exponentially leading asymptotic series near the real axis; 
it could not be trusted near the imaginary axis, as reqiU|red to get the spec- 
tral density. Quite recently, Dalmazi and Verbaarschot £3 have shown how to 
repair this latter deficiency of the large-mass expansion, by including other 
saddles. Indeed, their asymptotic expansion for large masses is valid also near 
the imaginary axis. Taking the discontinuity there according to eq. ( p^ ) pre- 
cisely reproduces the asymptotic expansion of the mircoscopic spectral density 
of the Dirac operator as computed based on the Random Matrix Theory for- 
mulation. This is a highly non-trivial fact, since in detail the computations 
are entirely different from those of either the Random Matrix Theory context 
or the supersymmetric formulation. Conceptually, the replica method has the 
advantage that the pattern of chiral symmetry breaking safely can be assumed 
to be the conventional one: S\J{Nf -1- iV„)L x SU(7V/ 4- Ny)R SU(7V/ -I- iV„), 
whereas in the supersymmstric formulation the symmetry breaking pattern is 
simply assumed to mimic as closely as possible the known bosonic one. The 
corresponding difficulty in the replica method is obviously how this pattern of 
chiral symmetry breaking can be given meaning for non-integer 7V„ (or rather, 
infinitesimal Ny , which is all that is required) . The fact that the replica method 
happens to agree with the supersymmetric one can, since they are so entirely 
different, be seen as an independent confirmation of both. 

4 Smallest Dirac operator eigenvalues near Tc 

Because of the Banks-Casher relation between the infinite-volume chiral con- 
densate S and the spectral density of the Dirac operator at the origin, E — 
7rp(0), it is obviously of interest to trace the depletion of Dirac eigenvalues at 
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aitc temperature. A lattice study of this question was first performed-in ref. 
, and this spring similar issues were addressed with higher statistics t3. 
The obvious question to ask is: can the Random Matrix Theory formu- 
lation of the effective Lagrangian (|l|) be used to predict the Dirac operator 
spectrum near Tc? The answer is clearly in the negative, since the effective 
finite-volume chiral Lagrangian (|l|) simply is incorrect at finite temperature. 
The best one can do is to take into account the leading effect of replacing the 
symmetric euclidean four- volume V = hy an asymmetric one, V = /T , 
where T still is on the order of 1/L. This means that this effective chiral La- 
grangian is suitable only for probing infinitesimally small temperatures, and 
the same goes for the associated Random Matrix Theory formulation. 

This does not mean that it is uninteresting to study the behavior of the 
smallest Dirac operator eigenvalues as the temperature is increased towards 
the critical temperalure Tc of the chiral phase transition. For instance, one 
observation of ref. E3 was that tracing the evolution in the magnitude of just 
the single smallest Dirac operator eigenvalue is a remarkably simple way to get 
an accurate detrmination of the phase transition point. Shown in figure 3 is 
an example of this, taken from ref. Ila. 

0.06 I • ^ • ^ • ^ • ^ • ^ 1 



0.04 



0.02 




Fig. 3 Just one single Dirac operator eigenvalue, traced through a series of ensembles with 
almost no averaging, suffices for getting an acciirate determination of the coupling at which 
the chiral phase transition occurs. Prom ref. t3. 
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A more ambituous goal would be to try to compute the corresponding 
microscopic spectral density at, say, precisely the critical temperature Tc- (For 
this to make sense, one should need a number of fermuions Nf for which the 
transition is continuous). The procedure is as clear as in the case of zero 
temperature, but a systematic way of obtaining the effective theory is lacking. 
Presumable the closest one can get is an effective Lagrangian framework like 
that of Pisarski and Wilczekcj. This concerns, however, an effective theory of 
which not even the leading behavior as a function of mass is known analytically. 
Progress on this front is more likely to come solely from the numerical side. 
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